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Neutron scattering experiments at the magnetic vacancy percolation threshold concentration, x v , 
using the random-field Ising crystal Feo.7eZno.24F2, show stability of the transition to long-range 
order up to fields H = 6.5 T. The observation of the stable long-range order corroborates the sharp 
boundary observed in computer simulations at x v separating equilibrium critical scattering behavior 
at high magnetic concentration from low concentration hysteretic behavior. Low temperature H > 
scattering line shapes exhibit the dependence on the scattering wavevector expected for percolation 
threshold fractal structures. 

PACS numbers: 



The dilute, anisotropic antiferromagnet Fe x Zni_ x F2 is 
an extensively studied prototype of the three dimensional 
(d = 3) random-field Ising model (RFIM) [§. As a result 
of the magnetic vacancies, the magnetic moment is not 
uniform and this allows a strong coupling to to an exter- 
nal magnetic field applied along the spin-ordering direc- 
tion. This constitutes the mechanism for the generation 
of random fields 0. It was shown that such a system 
is in the same universality class as a pure Ising magnet 
with random fields imposed 0. Settling the question of 
universality of the phase transition does not, however, 
address the effect of vacancies on microscopic domain 
formation, which can mask the phase transition in scat- 
tering experiments. Such microdomain formation, which 
occurs since domain walls can take advantage energet- 
ically of the vacancies, needs to be well understood in 
order to properly interpret the RFIM behavior of dilute 
magnets. 

For many years controversy surrounded the interpreta- 
tion of neutron scattering experiments jlj on the RFIM 
critical behavior of dilute anisotropic antiferromagnets 
in external magnetic fields, H, particularly Fe x Zni_ x F2 
and its less anisotropic isomorph Mn x Zni_ x F2. All of 
these studies, regardless of whether traditional scaling or 
various phcnomenological models were used in the inter- 
pretations, were done at concentrations x < 0.75 j|, |). 
This was natural since the strength of the random field 
increases with dilution and available field strengths re- 
quired high vacancy concentrations to readily create suit- 
ably strong random fields. It was, of course, realized that 
no ordering would take place for magnetic concentration 
below the magnetic percolation threshold concentration, 
x. p = 0.246. The magnetic vacancy percolation thresh- 
old concentration occurs at x v — 1 — x p = 0.754. Below 
this concentration, vacancies form a cluster that spans 
the crystal. The significant role of magnetic vacancy 
percolation in the formation of microdomains was not 



fully appreciated until recently |6j and, prior to that, it 
was widely assumed that microdomain formation was an 
intrinsic property of the RFIM as realized in dilute an- 
tiferromagnets. The microdomain structures for small 
x have been studied extensively ||. It has recently 
been shown that these structures play a crucial role in 
exchange-bias structures important to magnetic record- 
ing technology |j, Microscopic domain structure, for 
which the characteristic length scale is small compared to 
the instrumental resolution, masks the neutron scatter- 
ing critical behavior for two reasons. First, the scattering 
contribution from microscopic domains is superimposed 
on the scattering from thermal fluctuations, making it 
futile to separate the two. Second, there is a concomi- 
tant decrease in the Bragg scattering, which consequently 
no longer represents the strength of the RFIM order pa- 
rameter. This has been particularly frustrating, since 
characterization of the RFIM transition is important in 
light of the present disagreements between simulations 
and experiments ||. 

The critical behavior of Feo.93Zno.07F2 using neutron 
scattering techniques provided evidence [Tl| that micro- 
scopic domains could be avoided altogether by doing 
measurements at high magnetic concentrations, although 
very high quality crystals and high fields are required. 
Further experiments have been done using Feo.85Zno.15F2 
and Feo.87Zno.13F2 |f^, |f3). These experiments are pro- 
viding the avenue for a complete experimental character- 
ization of the RFIM universal critical behavior. It has 
become quite clear that the behavior at large x is quite 
distinct from that at low x which exhibits microdomain 
structure. Computer simulations || were done to model 
the behavior of the formation of microdomains and long- 
range order in Fc x Zni_ x F2 in an attempt to understand 
how the behavior crosses from one type of behavior to 
the other. It is suggested by these simulations that low- 
temperature metastability and microscopic domain for- 
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mation vanish abruptly above x = 0.76, which closely co- 
incides with the concentration of the magnetic vacancy 
percolation threshold concentration, x v = 0.754. Ap- 
parently, the percolating lattice of vacancies results in 
the instability of long-range order below the transition. 
In previous experiments, little attention has been paid 
to the percolation of magnetic vacancies. In light of its 
importance to the understanding of the RFIM we were 
motivated to investigate the scattering in Feo.75Zno.24F2, 
which is very close to the critical concentration x v . The 
concentration was determined using density measure- 
ments and the concentration gradient of a few tenths of 
a percent was determined using room temperature bire- 
fringence techniques [Q. 

Considerable focus has been given to the study of be- 
havior near the complementary threshold concentration 
for magnetic percolation, x p = 1 — x v , using neutron 
scattering, specific heat, linear birefringence, magnetiza- 
tion and ac susceptibility techniques |l5|, |l(| [lj], 18, p^|. 



The value x p — 0.246 is based on a calculation includ- 
ing only the dominant J2 interaction pTt . However, 
for concentrations close to x p , the system is extremely 
sensitive to very weak interactions that are insignifi- 
cant away from x p . Spin-glass- like behavior at H = 
has been well characterized. Even far above x p , large 
fields cause a crossover from the low-field microdomain- 
dominatcd random-field behavior to the spin-glass-like 
" 22l " 
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23] £4j, £5J £6[. 
The neutron scattering experiments were performed 
at the Oak Ridge National Laboratory High Flux Iso- 
tope Reactor using a double-axis spectrometer configu- 
ration. The beam was horizontally collimated to 20 min 
of arc before and after the sample and 48 min of arc be- 
fore the monochromator. The neutron energy was either 
13.7 meV on the HB1 spectrometer or 14.7 meV on the 
HB1A spectrometer. Higher energy neutrons were elim- 
inated using pyrolytic graphite filters. Most of the data 
were taken with transverse scans about the (100) antifer- 
romagnetic Bragg point. The Feo.76Zno.24F2 crystal has 
an irregular shape approximately 4 x 5 x 10 mm. It has a 
resolution limited Bragg peak, but very small secondary 
peaks appear for q > in the low temperature scans. 
Near the transition, these tiny peaks are not evident. All 
the data used in the analysis of the line shapes at low tem- 
peratures are on the q < side of the Bragg peak, where 
no hint of any secondary peaks are observed. The ther- 
mometry was based on a commercially calibrated carbon 
thermometer. Two primary thermal cycling procedures 
that are often employed to investigate hysteresis in the 
RFIM include: 1) cooling in the absence of a field, rais- 
ing the field and warming through the transitions (ZFC); 
and 2) cooling in the field (FC). 

Figure 1 shows scattering intensity vs. q, in reciprocal 
lattice units (rlu), at H = 3 and 5 T close to the tran- 
sition temperatures T c = 61.6 and 60.4 K, respectively. 
The transition at H = is at 63.2 K. Whereas the criti- 
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FIG. 1: ZFC and FC neutron scattering at H = 3 and 5 T 
below but close to the transition temperature. The contribu- 
tions include a constant background, critical scattering which 
is seen to increase as the temperature approaches (T C (H)), 
and a Bragg peak due to long range order. 



cal scattering from samples with x < x v exhibits strong 
hysteresis, it is clear that the |g| > critical scattering 
shown in Fig. 1 is free from hysteresis and that indicates 
there is no microscopic domain structure frozen in upon 
FC. 

For experiments free of extinction effects, the magnetic 
Bragg scattering intensity is expected to follow the power 
law behavior 
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where M s is the staggered magnetization and j3 « 0.35 
for the random-exchange model and (3 « 0.16 the RFIM 
Jl3| . However, neutron scattering in high-quality bulk 
crystals can suffer from severe extinction; the beam is de- 
pleted of neutrons that satisfy the Bragg condition and 
the scattering intensity is therefore saturated and can- 
not exhibit the correct T dependence. The extinction 
effects usually preclude determination of a reliable value 
of (3, the critical exponent for the staggered magneti- 
zation, from an analysis of the neutron scattering data 
in very high quality bulk crystals. The Bragg scatter- 
ing does show hysteresis, which indicates incomplete FC 
ordering on very long length scales, relative to the instru- 
mental resolution. Such hysteresis occurs for all x and 
is likely a consequence of the slow activated dynamics 
|7|, ||, |29|, |3(J of the RFIM very close to T C (H). 
The ZFC Bragg intensity, corrected for the background 
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FIG. 2: The ZFC antiferromagnetic Bragg peak for H = 3 T. 
The contribution at q = due to critical fluctuations is sub- 
tracted from the data represented by the open symbols to give 
the corrected intensities represented by the solid symbols. 



determined at high T, is shown in Fig. 2 vs. T for 
H = 3 T. The inset of Fig. 2 shows the critical scattering 
at small q, but well outside the transverse instrumental 
resolution where the Bragg intensity is negligible. The 
scattering line shapes in this sample are complicated by 
admixture of critical scattering and contributions from 
the vacancy lattice because of the proximity of this sam- 
ple to will be discussed. Hence, we could not 
confidently analyze in detail the critical scattering line 
shape. Nevertheless, we obtained an approximate ac- 
counting of the critical scattering intensity by taking a 
squared-Lorcntzian line shape folded with the appropri- 
ate resolution correction. An overall amplitude had to be 
chosen for the fit shown in the inset. This same amplitude 
was applied to the q = case and the resulting curve is 
shown in the main part of Fig. 2 as the solid curves. Sub- 
tracting this from the raw data (open symbols) yields the 
corrected Bragg scattering data (filled symbols). Taking 
into account the concentration gradient rounding of a few 
tenths of a percent, it is quite apparent that the Bragg 
intensity data approach T C (H) with a steep slope, in con- 
trast with scattering experiments |H) with x < x v , where 
the slope is nearly zero. Although we cannot analyze the 
data according to Eq. 1 to obtain (3 because of extinc- 
tion, we may still conclude from the shape of the Bragg 
intensity vs. T that this sample does not form micro- 
scopic domain structure. This is consistent with the lack 
of hysteresis in the critical scattering shown in Fig. 1. 



FIG. 3: The logarithm of the neutron scattering intensity vs. 
the logarithm of q for H = T and for ZFC and FC up to 
H = 7 T at T = 20K, well below T C (H). The background 
determined at large q has been subtracted. Only a few sets 
of data are shown for clarity. All ZFC data lie on the same 
curve, but only H = 5 and 7 T data are shown. 

At the percolation threshold concentrations, magnetic 
sites or magnetic vacancy sites form fractal structures. 
In either case, scattering from the fractal structure will 
exhibit a power law behavior (3^, |3|] 

■ (2) 

The only difference is that in the case of magnetic va- 
cancies there is also a Bragg scattering peak from the 
average M s . With magnetic site percolation, the aver- 
age M s is zero at the threshold. Since we believe the 
x = 0.76 crystal is close to x v , we plotted the logarithm 
of the scattering intensity vs. the logarithm of — q for 
q < in Fig. 3. Only some of the scans are shown; data 
were taken for H = and for 2T < H < 7 T in steps of 
0.5 T. Several interesting features are evident. 

As shown explicitly for two fields, H = 5 and 7 T in 
Fig. 3, the ZFC line shapes, for T = 29 K or less, are 
all identical with the H = line shape and are the low- 
est in intensity. For comparison, a line with a fractal 
exponent of 2.53 for three dimension, with the spectrom- 
eter resolution folded in, is plotted in the graph with the 
amplitude adjusted to fit the ZFC intensity data. It is 
clear that the ZFC scattering line shapes for x = 0.76 
follow Eq. 2 quite well. To contrast this behavior, we 
show in Fig. 3 similar data for a sample with concentra- 
tion of x w 0.87 indicated by the solid triangles, 
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for T = 58.5 K, only 3.4 K below the transition. The 
scattering for this concentration, well above x v , shows 
little evidence of scattering outside the Bragg region as 
expected since the vacancies do not form large fractal 
structures at this concentration. The behavior for ZFC 
shown in Fig. 3 suggests strongly that, for x — 0.76, the 
scattering is indeed from the vacancy percolation fractal 
structure under the ZFC procedure. 

The FC data for x = 0.76 increase in intensity with the 
applied field. We compare the data to Eq. 2 by adjust- 
ing the amplitude to fit the data at q = 10~ 1,7 rlu. It is 
quite clear that for H > the line shapes deviate strongly 
from the behavior in Eq. 2, more so as the field increases. 
The inset of Fig. 3 shows the deviations of the intensi- 
ties from the curves representing Eq. 2 at q — —0.007, 
—0.009, and —0.012 rlu as a function of the applied field. 
The deviations for H < 6.5 T increase smoothly with the 
field. Two possible sources exist for the excess scattering. 
One is the relief of extinction. This has been observed 
for neutron Bragg scattering in the RFIM experiments on 
bulk crystals (I] , but not for scattering outside the Bragg 
region. The other possibility, perhaps more significant, 
is the scattering from domains, which coexist with anti- 
ferromagnetic long-range order, that increase in number 
with increasing applied field. Since the total scattering 
in FC comes from sources in addition to that of vacancy 
sites, it is difficult to analyze in detail. 

The difference between FC intensities and Eq. 2 at 
H = 7 T deviate strongly from the smooth curves de- 
scribing the data for H < 6.5 T, as shown in the inset of 
Fig. 3. This qualitatively new behavior most likely rep- 
resents a breakdown of the antiferromagnetic long-range 
structure for H > 6.5 T. In such a case many more do- 
mains are introduced into the system, resulting in much 
more scattering intensity. This field is consistent with the 
increasing field at which spin-glass-like behavior appears 
for samples with x < x v , as shown in Fig. 4. Apparently, 
even at low temperatures, long-range antiferromagnetic 
order for x — 0.76 is stable upon FC for H < 7 T, in 
contrast to the behavior for x < x v , where metastable 
domains dominate the scattering under the FC proce- 
dure. For x w 0.87, only a few percent above Xy , mag- 
netization experiments fl34| indicate that the transition 
to long-range antiferromagnetic order breaks down only 
for fields above H — 18 T, nearly three times the field 
that causes a breakdown in the behavior of the sample 
with x = 0.76, demonstrating the stability of antiferro- 
magnetic order for x > x v . 

We have shown that the magnetic concentration range 
for which equilibrium random field critical scattering is 
observed for dilute antiferromagnets in an external field 
has a lower bound at x v . Whereas our x = 0.76 crystal 
shows no critical scattering hysteresis, it is quite evident 
for slightly smaller magnetic concentrations ||. Our re- 
sults suggest that the percolation of vacancy sites occur- 
ring for x < x v precipitates the formation of domains 
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FIG. 4: The concentration dependence of the phases observed 
in Fe x Zni_ x F2. For x < x p , the magnetic percolation thresh- 
old concentration, only antiferromagnetic short-range order 
is possible. For x p < x < x v , microscopic domain forma- 
tion occurs at low fields and spin-glass-like behavior occurs 
at high fields. For x > x v , antiferromagnetic long-range or- 
der is stable, without microscopic domain formation, to very 
large applied fields. Critical behavior measurements can only 
be reliably done for x > x v . Data are taken from various 
experiments cited in the text. 

below T C (H) in ZFC preparation as well as in FC, cor- 
roborating the conclusions drawn from simulations Q . A 
theoretical connection between the one dimensional frac- 
tal geometry of vacancy sites at percolation and three 
dimensional domains has not been adequately explained 
from a theoretical perspective. 

It is now evident that there are three magnetic con- 
centration regimes in d = 3 dilute antiferromagnets sepa- 
rated by the percolation threshold concentrations x p and 
x v . For x > x v , long-range antiferromagnetic order is 
stable up to very large magnetic fields. For x p < x < x v 
the system is unstable. At low fields, the formation 
of microdomain structure takes place upon FC for all 
T < T C (H) and upon ZFC close to T C (H), A spin-glass- 
like phase forms at high field. Below x p , there can be 
no long-range magnetic order. It is the geometry of the 
lattice in question which defines the location of these 
boundaries, and although we study one particular mag- 
netic lattice type, the body centered tetragonal structure 
of Fe x Zni_ x F2, our results should apply more generally 
to dilute magnets in an applied field. 

Interestingly, the specific heat behavior is not depen- 
dent in an obvious way on the concentration. Similar 
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hysteresis upon FC and ZFC is observed [0 very close 
to T C (H) for concentrations above and below x v . No 
specific heat hysteresis is observed at low T. The con- 
trast between the relative insensitivity of the specific heat 
techniques with the extreme sensitivity of the scattering 
techniques is certainly due to the greater dependence of 
the scattering on long length correlations that are greatly 
affected by domain formation. The hysteresis in the case 
of specific heat is related to the activated dynamics very 
close to T C (H) that affects the behavior at all x and not 
domain formation, which only occurs for x < x v . 

From the results of this investigation, we conclude that 
studies of the random-field phase transition should be 
conducted with magnetic concentrations greater than x v . 
It is advantageous to use concentrations not too much 
greater than this to maximize the random-fields for avail- 
able applied fields. However, if the concentration is too 
close to x v , one must take into account scattering from 
the magnetic vacancy percolation cluster. Recent exper- 
iments at x — 0.87 indicate that, at this concentration, 
such scattering is negligible [|l2| [l| . 
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